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FIBER RENDERING APPARATUS 

BACKGROUND OF THE INVENTION 

The present invention relates to a fiber rendering method and MRI 
5 (magnetic resonance imaging) apparatus, and more particularly to a method and 
MRI apparatus for properly ^rendering brain white matter fibers obtained by 
diffusion tensor imaging. 

Figure 21 is a flow chart showing a conventional fiber rendering method. 

At Step PI, an MR image in an axial or oblique plane is produced from 
10 three-dimensional image data collected by a diffusion tensor method or another 
imaging method (Tl- or T2-enhanced or the like) in an MRI apparatus, and the 
MR image is displayed. 

At Step P2, an operator specifies a two-dimensional region of interest Rl 
(or a three-dimensional volumetric region of interest) on a displayed MR image 
1 5 Gl, as shown in Figure 22. 

At Step P3', regular grid points are generated in the region of interest Rl (or 
in the volumetric region of interest) as shown in Figure 23, and they are defined 
as tracking start points SI, S2, S3, — 

At Step P5, one of the tracking start points is selected. 
20 At Step P6', diffusion tensor analysis is performed on the selected tracking 

start point in the three-dimensional image data collected by the diffusion tensor 
method in the MRI apparatus to determine the direction of the principal axis 
vector, i.e., the direction of the first eigenvector. 

At Step P7, if a point at a unit distance along the direction of the principal 
25 axis vector falls within the three-dimensional image data space, the point is 
defined as a neighbor point and the flow proceeds to Step P8'; and if the point 
falls outside the three-dimensional image data space, the flow proceeds to Step 
Pll. 

At Step P8', data at the neighbor point is created by interpolation or the like 
30 on the three-dimensional image data, and diffusion tensor analysis is performed 
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to determine the direction of the principal axis vector and the FA (fractional 
anisotropy) value. 

At Step P9, if the FA value is equal to or more than a threshold, the flow 
goes back to Step P7 to continue the fiber tracking because the fiber tracking has 
5 not reached an end portion of a brain white matter fiber; and if the FA value is 
less than the threshold, the flow proceeds to Step Pll to terminate the fiber 
tracking because an end portion of a brain white matter fiber has been reached. 

In this way, Steps P7 — P9 are repeated until no more three-dimensional 
image data are found or the fiber tracking has reached an end portion of a brain 
10 white matter fiber, and a fiber is tracked from the tracking start point SI to a 
neighbor point Nl, N2, N3, . . ., as exemplarily shown in Figure 24. At that time, 
connectivity is decided by using a scalar product of vectors, for example. 

At Step Pll, points from the tracking start point to the last neighbor point 
are saved as one brain white matter fiber. 
15 At Step P12, if any tracking start point not selected at Step P5 remains, the 

flow goes back to Step P5; otherwise, proceeds to Step P14'. 

At Step P14', an image of the saved brain white matter fibers as viewed in a 
desired view direction is produced and displayed, as exemplarily shown in 
Figure 25. 

20 A diffusion tensor and a nerve fiber extending direction are described in, 

for example, "Microstructural and Physiological Features of Tissues Elucidated 
by Quantitative-Diffusion-Tensor MRI" by Peter J. Basser and Carlo Pierpaoli, 
Journal of Magnetic Resonance, Series B 111, pp. 209 — 219 (1996), and in 
"Dif fusion Anisotropy — 2D and 3D images of Brain White Matter Fibers — " by 

25 Yasuomi Kinosada (Kyoto Prefectural University of Medicine, Department of 
Radiology), the 30 th Meeting of MR Imaging Study Group, September 4, 1998, at 
Sapporo, Japan. 

SUMMARY OF THE INVENTION 
30 When the grid points regularly generated at Step P3 in Figure 21 are 



defined as the tracking start points, the fiber density looks as if it suddenly 
decreases when the view direction is parallel to a direction of the grid point 
arrangement, because the nerve fibers passing through the tracking start points 
lining up in the view direction appear to overlap one another* leading to a 
problem that the image gives an unnatural impression. 

Thus, a first object of the present invention is to provide a fiber rendering 
method capable of preventing a situation in which the fiber density looks as if it 
suddenly decreases in a specific view direction. 

When the threshold at Step P9 in Figure 21 is small, even a portion having a 
considerably low FA value, i.e., a portion with considerably low fiber tracking 
reliability, is rendered. The portion with considerably low fiber tracking 
reliability is, however, rendered in the same manner of display as that for 
rendering a portion with high reliability, and these portions cannot be 
distinguished, leading to the problem that this poses a hindrance to accurate 
diagnosis. On the other hand, when the threshold at Step P9 in Figure 21 is 
large, the fiber tracking is aborted before an end portion of a brain white matter 
fiber is reached, leading to a problem that fibers cannot be fully rendered. 

Thus, a second object of the present invention is to provide a fiber 
rendering method capable of rendering fibers in a manner of display that 
incorporates the degree of fiber tracking reliability. 

In the conventional technique, since eigenvalues of diffusion tensors are 
not incorporated in display of tracked fibers, there is a problem that variation in 
eigenvalues of diffusion tensors cannot be seen when observing the rendered 
fibers. 

Thus, a third object of the present invention is to provide a fiber rendering 
method capable of rendering fibers in a manner of display that incorporates 
variation in the eigenvalues of diffusion tensors. 

As shown in Figure 27, at a nerve fiber intersection C, nerve fibers having 
different connection directions intersect each other. The conventional tracking, 
however, employs only the direction of the principal axis vector at a selected 
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neighbor point, and therefore it cannot distinguish between fibers that intersect 
each other at a nerve fiber intersection, leading to a problem that the tracking 
direction may be mistaken as shown in Figure 28. 

Thus, a fourth object of the present invention is to provide a fiber rendering 
5 method capable of conducting tracking without mistaking the direction even at a 
portion where nerve fibers having different connection directions intersect each 
other. 

In diagnosing leukodystrophy, for example, knowledge about whether 
connection by fiber nerves between two sites has been destroyed provides useful 
10 information. 

Thus, a fifth object of the present invention is to provide a fiber rendering 
method by which connectivity by fiber nerves between two sites that an operator 
specifies can be visually recognized. 

In accordance with its first aspect, the present invention provides a fiber 

15 rendering method characterized in comprising: specifying a region of interest or 
volumetric region of interest in three-dimensional image data collected by a 
diffusion tensor method in an MRI apparatus; defining regular grid points in the 
region of interest or volumetric region of interest; then defining points obtained 
by randomly moving the grid points in a two-dimensional or three-dimensional 

20 manner as tracking start points; performing diffusion tensor analysis on each 
tracking start point in the three-dimensional image data to determine a direction 
of a principal axis vector; tracking a fiber by repeatedly selecting a neighbor 
point along the direction of the principal axis vector and performing diffusion 
tensor analysis on the neighbor point to determine a direction of a principal axis 

25 vector; and producing and displaying an image of the tracked fibers as viewed in 
a desired view direction. 

In the fiber rendering method of the first aspect, the number of tracking 
start points overlapping one another is approximately the same in any view 
direction. Therefore, a situation in which the fiber density looks as if it 

30 suddenly decreases in a specific view direction is prevented. Taking an overall 
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view of the region of interest or volumetric region of interest, the density of the 
track start points is uniform and no density variation occurs. 

In accordance with its second aspect, the present invention provides a fiber 
rendering method characterized in comprising: defining tracking start points in 
5 three-dimensional image data collected by a diffusion tensor method in an MRI 
apparatus; performing diffusion tensor analysis on each tracking start point in 
the three-dimensional image data to determine a direction of a principal axis 
vector and a diffusion anisotropy value; tracking a fiber by repeatedly selecting a 
neighbor point along the direction of the principal axis vector and performing 

10 diffusion tensor analysis on the neighbor point to determine a direction of a 
principal axis vector and a diffusion anisotropy value; and producing an image 
of the tracked fibers as viewed in a desired view direction and displaying the 
image with opacity reflecting the diffusion anisotropy values at the tracking start 
points and neighbor points. 

15 In the fiber rendering method of the second aspect, the transparency of a 

fiber to be rendered is changed according to the diffusion anisotropy value. 
Therefore, the degree of fiber tracking reliability can be visually recognized from 
the transparency of the rendered fibers. 

In accordance with its third aspect, the present invention provides the fiber 

20 rendering method having the aforementioned configuration, characterized in 
that an FA value is used as the diffusion anisotropy value. 

In the fiber rendering method of the third aspect, the transparency of a fiber 
to be rendered can be changed according to an FA value that takes a value 
between zero and one depending upon the diffusion anisotropy. 

25 In accordance with its fourth aspect, the present invention provides the 

fiber rendering method having the aforementioned configuration, characterized 
in that: 

Xn+l ~ FA n * X n/ 

where X n +i represents an opacity at a neighbor point, FA n represents an FA value 
30 at the immediately preceding neighbor point or tracking start point, and X n 
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represents an opacity thereat. 

In the fiber rendering method of the fourth aspect, the transparency can be 
gradually increased from the tracking start point toward an end portion, and 
sharply increased at the end portion. 
5 In accordance with its fifth aspect, the present invention provides a fiber 

rendering method characterized in comprising: defining tracking start points in 
three-dimensional image data collected by a diffusion tensor method in an MRI 
apparatus; performing diffusion tensor analysis on each tracking start point in 
the three-dimensional image data to determine a direction of a principal axis 
10 vector and eigenvalues of a diffusion tensor; tracking a fiber by repeatedly 
selecting a neighbor point along the direction of the principal axis vector and 
performing diffusion tensor analysis on the neighbor point to determine a 
direction of a principal axis vector and eigenvalues of a diffusion tensor; and 
producing an image of the tracked fibers as viewed in a desired view direction 
15 and displaying the image with display colors reflecting the eigenvalues of the 
diffusion tensors at the tracking start points and neighbor points. 

In the fiber rendering method of the fifth aspect, the display color of fibers 
to be rendered is changed according to the eigenvalues of the diffusion tensors. 
Therefore, the change in the eigenvalues of the diffusion tensors can be visually 
20 recognized by the change in the display color of the rendered fibers. 

In accordance with its sixth aspect, the present invention provides the fiber 
rendering method having the aforementioned configuration, characterized in 
that: a display color (R, G, B) is defined as: 
R:G:B = 1:X2/M:A3 / XI, 
25 where XI, X2 and A3 represent eigenvalues of a diffusion tensor. 

In the fiber rendering method of the sixth aspect, the diffusion can be 
known to be more isotropic as the display color is closer to white, and to be more 
anisotropic as the display color is closer to red. 

In accordance with its seventh aspect, the present invention provides a 
30 fiber rendering method characterized in comprising: defining tracking start 
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points in three-dimensional image data collected by a diffusion tensor method in 
an MRI apparatus; performing diffusion tensor analysis on each tracking start 
point in the three-dimensional image data to determine a direction of a principal 
axis vector and defining the direction of the principal axis vector as a tracking 
5 direction vector; tracking a fiber by repeatedly selecting a neighbor point along 
the tracking direction vector, performing diffusion tensor analysis on the 
neighbor point to obtain diffusion tensor information, and deterirrining a 
tracking direction vector from the diffusion tensor information and at least an 
immediately preceding tracking direction vector; and producing and displaying 

10 an image of the tracked fibers as viewed in a desired view direction. 

In the fiber rendering method of the seventh aspect, since a new tracking 
direction vector is determined from diffusion tensor information of a neighbor 
point and at least an immediately preceding tracking direction vector, nerve 
fibers in different connection directions can be distinguished based on the 

15 preceding connection directions even at a portion at which the nerve fibers in 
different connection directions intersect each other, and the nerve fibers can be 
tracked without mistaking the direction. 

In accordance with its eighth aspect, the present invention provides the 
fiber rendering method having the aforementioned configuration, characterized 

20 in that: 

d i+ i = {XI (ei - di) ei + X2 (e 2 • di) e 2 + A3 (e 3 • di) e 3 } 
/ | XI (ei ■ di) ei + X2 (e 2 • di) e 2 + X3 (e 3 • di) e 3 | , 
where XI, X2 and A3 represent eigenvalues of a diffusion tensor at a neighbor 
point, el, e2 and e3 represent eigenvectors thereat, di+i represents a tracking 
25 direction vector thereat, and di represents a tracking direction vector at an 
immediately preceding neighbor point or tracking start point. 

In the fiber rendering method of the eighth aspect, a tracking direction 
vector di+i can be determined from an immediately preceding tracking direction 
vector di, and eigenvalues of a diffusion tensor Al, X2 and A3 and eigenvectors el, 
30 e2 and e3 at a neighbor point. 
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In accordance with its ninth aspect, the present invention provides a fiber 
rendering method characterized in comprising: specifying a start region of 
interest and an end region of interest or a start volumetric region of interest and 
an end volumetric region of interest in three-dimensional image data collected by 
5 a diffusion tensor method in an MRI apparatus; defining tracking start points in 
the start region of interest or start volumetric region of interest; tracking a fiber 
by performing diffusion tensor analysis from each tracking start point in the 
three-dimensional image data; deciding whether each tracked fiber passes 
through the end region of interest or end volumetric region of interest; and 

10 producing and displaying an image of only the fibers that are decided to pass 
through as viewed in a desired view direction. 

In the fiber rendering method of the ninth aspect, since only the nerve 
fibers passing through two sites are rendered, connectivity of the nerve fibers 
between the two sites can be visually recognized. 

15 In accordance with its tenth aspect, the present invention provides the fiber 

rendering method having the aforementioned configuration, characterized in 
comprising: calculating and displaying a total sum with respect to all the fibers 
decided to pass through: 

M_Value = S XI - FA / L, 

20 where XI represents a first eigenvalue of a diffusion tensor of a fiber decided to 
pass through, FA represents an FA value thereof, and L represents the total 
length of the fiber. 

In the fiber rendering method of the tenth aspect, quantitative assessment is 
enabled by employing M_Value as an indicator of the strength of connection by 
25 nerve fibers between two sites. 

In accordance with its eleventh aspect, the present invention provides a 
fiber rendering apparatus characterized in comprising: means for specifying a 
region of interest or volumetric region of interest in three-dimensional image 
data collected by a diffusion tensor method in an MRI apparatus; means for 
30 defining regular grid points in the region of interest or volumetric region of 
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interest; means for defining points obtained by randomly moving the grid points 
in a two-dimensional or three-dimensional manner as tracking start points; 
means for performing diffusion tensor analysis on each tracking start point in the 
three-dimensional image data to determine a direction of a principal axis vector; 
5 means for tracking a fiber by repeatedly selecting a neighbor point along the 
direction of the principal axis vector and performing diffusion tensor analysis on 
the neighbor point to determine a direction of a principal axis vector; and means 
for producing and displaying an image of the tracked fibers as viewed in a 
desired view direction. 

10 In the fiber rendering apparatus of the eleventh aspect, the fiber rendering 

method of the first aspect can be suitably implemented. 

In accordance with its twelfth aspect, the present invention provides a fiber 
rendering apparatus characterized in comprising: means for defining tracking 
start points in three-dimensional image data collected by a diffusion tensor 

15 method in an MRI apparatus; means for performing diffusion tensor analysis on 
each tracking start point in the three-dimensional image data to determine a 
direction of a principal axis vector and a diffusion anisotropy value; means for 
tracking a fiber by repeatedly selecting a neighbor point along the direction of the 
principal axis vector and performing diffusion tensor analysis on the neighbor 

20 point to determine a direction of a principal axis vector and a diffusion 
anisotropy value; and means for producing an image of the tracked fibers as 
viewed in a desired view direction and displaying the image with opacity 
reflecting the diffusion anisotropy values at the tracking start points and 
neighbor points. 

25 In the fiber rendering apparatus of the twelfth aspect, the fiber rendering 

method of the second aspect can be suitably implemented. 

In accordance with its thirteenth aspect, the present invention provides the 
fiber rendering apparatus having the aforementioned configuration, 
characterized in that an FA value is used as the diffusion anisotropy value. 

30 In the fiber rendering apparatus of the thirteenth aspect, the fiber rendering 
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method of the third aspect can be suitably implemented. 

In accordance with its fourteenth aspect, the present invention provides the 
fiber rendering apparatus having the aforementioned configuration, 
characterized in that: 

5 X n +1 — FA n • Xn, 

where X n +i represents an opacity at a neighbor point, FA n represents an FA value 
at the immediately preceding neighbor point or tracking start point, and X n 
represents an opacity thereat. 

In the fiber rendering apparatus of the fourteenth aspect, the fiber 

10 rendering method of the fourth aspect can be suitably implemented. 

In accordance with its fifteenth aspect, the present invention provides a 
fiber rendering apparatus characterized in comprising: means for defining 
tracking start points in three-dimensional image data collected by a diffusion 
tensor method in an MRI apparatus; means for performing diffusion tensor 

15 analysis on each tracking start point in the three-dimensional image data to 
determine a direction of a principal axis vector and eigenvalues of a diffusion 
tensor; means for tracking a fiber by repeatedly selecting a neighbor point along 
the direction of the principal axis vector and performing diffusion tensor analysis 
on the neighbor point to determine a direction of a principal axis vector and 

20 eigenvalues of a diffusion tensor; and means for producing an image of the 
tracked fibers as viewed in a desired view direction and displaying the image 
with display colors reflecting the eigenvalues of the diffusion tensors at the 
tracking start points and neighbor points. 

In the fiber rendering apparatus of the fifteenth aspect, the fiber rendering 

25 method of the fifth aspect can be suitably implemented. 

In accordance with its sixteenth aspect, the present invention provides the 
fiber rendering apparatus having the aforementioned configuration, 
characterized in that: a display color (R, G, B) is defined as: 
R:G:B = 1:A2/A1:A3/ Al, 

30 where Al, A2 and A3 represent eigenvalues of a diffusion tensor. 
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In the fiber rendering apparatus of the sixteenth aspect, the fiber rendering 
method of the sixth aspect can be suitably implemented. 

In accordance with its seventeenth aspect, the present invention provides a 
fiber rendering apparatus characterized in comprising: means for defining 
5 tracking start points in three-dimensional image data collected by a diffusion 
tensor method in an MRI apparatus; means for performing diffusion tensor 
analysis on each tracking start point in the three-dimensional image data to 
determine a direction of a principal axis vector and defining the direction of the 
principal axis vector as a tracking direction vector; means for tracking a fiber by 
10 repeatedly selecting a neighbor point along the tracking direction vector, 
performing diffusion tensor analysis on the neighbor point to obtain diffusion 
tensor information, and determining a tracking direction vector from the 
diffusion tensor information and at least an immediately preceding tracking 
direction vector; and means for producing and displaying an image of the 
15 tracked fibers as viewed in a desired view direction. 

In the fiber rendering apparatus of the seventeenth aspect, the fiber 
rendering method of the seventh aspect can be suitably implemented. 

In accordance with its eighteenth aspect, the present invention provides the 
fiber rendering apparatus having the aforementioned configuration, 
20 characterized in that: 

d i+1 = {XI (ei - di) ei + X2 (e 2 • di) e 2 + 73 (e 3 - di) e 3 } 
/ | XI (ei - di) ei + X2 (e 2 - di) e 2 + A3 (e 3 • di) e 3 1 , 
where XI, X2 and A3 represent eigenvalues of a diffusion tensor at a neighbor 
point, el, e2 and e3 represent eigenvectors thereat, di+i represents a tracking 
25 direction vector thereat, and di represents a tracking direction vector at an 
immediately preceding neighbor point or tracking start point. 

In the fiber rendering apparatus of the eighteenth aspect, the fiber 
rendering method of the eighth aspect can be suitably implemented. 

In accordance with its nineteenth aspect, the present invention provides a 
30 fiber rendering apparatus characterized in comprising: means for specifying a 
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start region of interest and an end region of interest or a start volumetric region 
of interest and an end volumetric region of interest in three-dimensional image 
data collected by a diffusion tensor method in an MRI apparatus; means for 
defining tracking start points in the start region of interest or start volumetric 
5 region of interest; means for tracking a fiber by performing diffusion tensor 
analysis from each tracking start point in the three-dimensional image data; 
means for deciding whether each tracked fiber passes through the end region of 
interest or end volumetric region of interest; and means for producing and 
displaying an image of only the fibers that are decided to pass through as viewed 

10 in a desired view direction. 

In the fiber rendering apparatus of the nineteenth aspect, the fiber 
rendering method of the ninth aspect can be suitably implemented. 

In accordance with its twentieth aspect, the present invention provides the 
fiber rendering apparatus having the aforementioned configuration, 

15 characterized in comprising: means for calculating and displaying a total sum 
with respect to all the fibers decided to intersect: 

M_Value = E XI - FA / L, 
where XI represents a first eigenvalue of a diffusion tensor of a fiber decided to 
pass through, FA represents an FA value thereof, and L represents the total 

20 length of the fiber. 

In the fiber rendering apparatus of the twentieth aspect, the fiber rendering 
method of the tenth aspect can be suitably implemented. 

According to the fiber rendering method and fiber rendering apparatus of 
the present invention, the following effects can be obtained: 

25 (1) The number of tracking start points overlapping one another is 
approximately the same in various view directions. Therefore, a situation in 
which the fiber density looks as if it suddenly decreases in a specific view 
direction is prevented. Taking an overall view of the region of interest or 
volumetric region of interest, the density of the track start points is uniform and 

30 no density variation occurs; 
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(2) A portion of a rendered fiber having low transparency can be known to 
have high fiber tracking reliability, and a portion having high transparency can 
be known to have low fiber tracking reliability. Therefore, even when a portion 
having considerably low fiber tracking reliability is rendered, the portion with 

5 considerably low fiber tracking reliability and the portion with high reliability 
can be distinguished, which avoids hindrance to accurate diagnosis; 

(3) Whether diffusion is isotropic or anisotropic can be visually recognized 
from the display color of rendered fibers; 

(4) Nerve fibers in different connection directions can be distinguished based 
10 on the preceding connection directions even at a portion at which the nerve fibers 

in different connection directions intersect each other, and the nerve fibers can be 
tracked without mistaking the direction; 

(5) Since only the nerve fibers passing through two sites can be rendered, 
connectivity of the nerve fibers between the two sites can be visually recognized; 

15 and 

(6) Quantitative assessment on the strength of connection by nerve fibers 
between two sites is enabled. 

Further objects and advantages of the present invention will be apparent 
from the following description of the preferred embodiments of the invention as 
20 illustrated in the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram showing an MRI apparatus in accordance with 
a first embodiment. 

25 Figure 2 is a flow chart showing fiber rendering processing in accordance 

with the first embodiment. 

Figure 3 is a flow chart continued from Figure 2. 

Figure 4 exemplarily shows a screen for specifying a region of interest. 
Figure 5 exemplarily shows regularly arranged grid points. 
30 Figure 6 exemplarily shows irregularly position-shifted tracking start 
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points. 

Figure 7 is a conceptual diagram showing a fiber tracking condition. 

Figure 8 exemplarily shows an image of obtained fibers as viewed in a 
desired view direction. 

Figure 9 exemplarily shows an image of the obtained fibers as viewed in 
another view direction. 

Figure 10 is a flow chart showing fiber rendering processing in accordance 
with a second embodiment. 

Figure 11 is a flow chart continued from Figure 10. 

Figure 12 is a flow chart continued from Figure 11. 

Figure 13 is a conceptual diagram showing a tracking direction vector. 

Figure 14 is a conceptual diagram showing that the tracking directions are 
not mistaken even if fibers intersect each other. 

Figure 15 is an explanatory diagram showing that the tracking direction is 
not mistaken even at a nerve fiber intersection. 

Figure 16 is a flow chart showing fiber rendering processing in accordance 
with a third embodiment. 

Figure 17 is a flow chart continued from Figure 16. 

Figure 18 is a flow chart continued from Figure 17. 

Figure 19 exemplarily shows a screen for specifying start and end regions 
of interest. 

Figure 20 exemplarily shows a screen that displays only the fibers 
connecting the start and end regions of interest. 

Figure 21 is a flow chart showing conventional fiber rendering processing. 

Figure 22 exemplarily shows a screen for specifying a region of interest. 

Figure 23 exemplarily shows regularly arranged tracking start points. 

Figure 24 is a conceptual diagram showing a fiber tracking condition. 

Figure 25 exemplarily shows an image of obtained fibers as viewed in a 
desired view direction. 

Figure 26 exemplarily shows an image of the obtained fibers as viewed in 



16NM02038 



-15- 

another view direction. 

Figure 27 is a conceptual diagram showing that fibers intersect each other 
at a nerve fiber intersection. 

Figure 28 is an explanatory diagram showing that the tracking direction is 
5 mistaken at the nerve fiber intersection. 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention will now be described in more detail with reference 
to embodiments shown in the accompanying drawings. 
10 - First Embodiment - 

Figure 1 is a block diagram showing an MRI apparatus in accordance with 
one embodiment of the present invention. 

In the MRI apparatus 100, a magnet assembly 1 has a bore (cavity portion) 
for inserting therein a subject, and is provided, surrounding the bore, with a 
15 gradient coil (which comprises X-axis, Y-axis and Z-axis coils, and the 
combination thereof determines slice, warp and read axes) 1G for generating 
gradient magnetic fields, a transmit coil IT for applying RF pulses for exciting 
spins of atomic nuclei within the subject, a receive coil 1R for detecting NMR 
signals from the subject, and a static magnetic field power supply 2 and static 
20 magnetic field coil 1C for generating a static magnetic field. 

It should be noted that permanent magnets may be employed in place of 
the static magnetic field power supply 2 and static magnetic field coil 1C 
(superconductive coil). 

The gradient coil 1G is connected to a gradient coil driving circuit 3. The 
25 transmit coil IT is connected to an RF power amplifier 4. The receive coil 1R is 
connected to a preamplifier 5. 

A sequence memory circuit 8 operates the gradient coil driving circuit 3 
based on a stored pulse sequence in response to instructions from a computer 7 
to thereby generate gradient magnetic fields from the gradient coil 1G. The 
30 sequence memory circuit 8 also operates a gate modulation circuit 9 to modulate 
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high frequency output signals from an RF oscillation circuit 10 into pulsed 
signals of predefined timing and envelope. The pulsed signals are applied to 
the RF power amplifier 4 as excitation pulses, power-amplified in the RF power 
amplifier 4, and then applied to the transmit coil IT in the magnet assembly 1 to 
5 transmit RF pulses. 

The preamplifier 5 amplifies NMR signals from the subject detected at the 
receive coil 1R in the magnet assembly 1, and inputs the signals to a phase 
detector 12. The phase detector 12 phase-detects the NMR signals from the 
preamplifier 5 employing the output from the RF oscillation circuit 10 as a 

10 reference signal, and supplies the phase-detected signals to an A/D converter 11. 
The A/D converter 11 converts the phase-detected analog signals into MR data in 
the form of digital signals, and inputs them to the computer 7. 

The computer 7 reads the MR data from the A/D converter 11, and 
performs image reconstruction calculation to produce an MR image. The 

15 computer 7 is also responsible for overall control such as receiving information 
supplied from an operator console 13. Furthermore, the computer 7 conducts 
fiber rendering processing, which will be described later with reference to Figure 
2. 

A display device 6 displays the MR image and a fiber image which will be 
20 described later. 

Figure 2 is a flow chart showing fiber rendering processing by the MRI 
apparatus 100. 

At Step PI, an MR image in an axial or oblique plane is produced from 
three-dimensional image data collected by a diffusion tensor method or another 
25 imaging method (Tl- or T2-enhanced or the like) in the MRI apparatus 100, and 
the MR image is displayed. 

At Step P2, an operator specifies a two-dimensional region of interest Rl 
(or a three-dimensional volumetric region of interest) on a displayed MR image 
Gl, as shown in Figure 4. 
30 At Step P3, regular grid points gl, g2, g3, ... are generated in the region of 
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interest Rl (or in the volumetric region of interest), as shown in Figure 5. 

At Step P4, points obtained by randomly moving the grid points gl, g2, g3, 
... in a two-dimensional (or three-dimensional) manner are defined as tracking 
start points SI, S2, S3, as shown in Figure 6. Random numbers for the 
5 random moving can be generated using a distribution function such as a 
Gaussian distribution or uniform distribution. The range of the moving may be 
defined so that most of the points after the moving fall within intervals between 
the grid points gl, g2, g3, . . . . 

At Step P5, one of the tracking start points is selected. 
10 At Step P6, diffusion tensor analysis is performed on the selected tracking 

start point in the three-dimensional image data collected by the diffusion tensor 
method in the MRI apparatus 100 to determine the direction of the principal axis 
vector, the FA value, and the eigenvalues. 

At Step P7, if a point at a unit distance along the direction of the principal 
15 axis vector falls within the three-dimensional image data space, the point is 
defined as a neighbor point and the flow proceeds to Step P8; and if the point 
falls outside the three-dimensional image data space, the flow proceeds to Step 
Pll. 

At Step P8, data at the neighbor point is created by interpolation or the like 
20 on the three-dimensional image data, and diffusion tensor analysis is performed 
to determine the direction of the principal axis vector, the FA value, and the 
eigenvalues. 

At Step P9, if the FA value is equal to or more than a threshold, the flow 
goes back to Step P7 to continue the fiber tracking because the fiber tracking has 
25 not reached an end portion of a brain white matter fiber; and if the FA value is 
less than the threshold, the flow proceeds to Step Pll to terminate the fiber 
tracking because an end portion of a brain white matter fiber has been reached. 

In this way, Steps P7 — P9 are repeated until no more three-dimensional 
image data are found or the fiber tracking has reached an end portion of a brain 
30 white matter fiber, and a fiber is tracked from the tracking start point SI to a 
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neighbor point Nl, N2, N3, as exemplarily shown in Figure 7. At that time, 
connectivity is decided by using a scalar product of vectors, for example. 

At Step Pll, points from the tracking start point to the last neighbor point 
are saved as one brain white matter fiber. 
5 At Step P12, if any tracking start point not selected at Step P5 remains, the 

flow goes back to Step P5; otherwise, proceeds to Step P14 in Figure 3. 

At Step P14 in Figure 3, an image of the saved brain white matter fibers as 
viewed in a desired view direction is produced, as exemplarily shown in Figure 
8. 

10 At Step P15, the opacity at the tracking start point is defined as Xo. 

Moreover, 

X n +1 = FA n • X n 

is set, where X n +i represents the opacity at a neighbor point, FA n represents the 
FA value at the immediately preceding neighbor point or tracking start point, 
1 5 and X n represents the opacity thereat. 

At Step PI 6, the display color (R, G, B) is defined as: 
R : G : B = 1 : X2 / XI : 73 / XI, 
where XI, X2 and A3 represent the eigenvalues of the diffusion tensor. 

At Step PI 7, the image of the fibers is displayed using the opacity X and the 
20 display color (R, G, B). 

According to the MRI apparatus 100 of the first embodiment, the following 
effects can be obtained: 

(1) As shown in Figures 8 and 9, the number of tracking start points 
overlapping one another is approximately the same in various view directions. 

25 Therefore, a situation in which the fiber density looks as if it suddenly decreases 
in a specific view direction is prevented. Taking an overall view of the region of 
interest or volumetric region of interest, the density of the track start points is 
uniform and no density variation occurs; 

(2) A portion of a rendered fiber having low transparency can be known to 
30 have high fiber tracking reliability, and a portion having high transparency can 
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be known to have low fiber tracking reliability. Therefore, even when a portion 
having considerably low fiber tracking reliability is rendered by decreasing the 
threshold at Step P9 in Figure 2, the portion with considerably low fiber tracking 
reliability and the portion with high reliability can be distinguished, which 
avoids hindrance to accurate diagnosis; and 

(3) The diffusion can be known as being more isotropic as the display color for 
the rendered fibers is closer to white, and as being more anisotropic as the 
display color is closer to red. 

In addition, modifications as follows may be made: 

(1) The opacity X may be calculated based on another indicator that reflects the 
diffusion anisotropy (for example, the eigenvalue ratio, A2 / XI, A3 / XI, relative 
anisotropy, volume ratio); and 

(2) The display color (R, G, B) may be determined as R : G : B = XI / (XI + X2 + 
13) : X2 / (XI + X2 + X3) : A3 / (A.1 + X2 + A3). 

- Second Embodiment - 

Figure 10 is a flpw chart showing fiber rendering processing by an MRI 
apparatus in accordance with a second embodiment. 

At Step Ql, an MR image in an axial or oblique plane is produced from 
three-dimensional image data collected by a diffusion tensor method or another 
imaging method (Tl- or T2-enhanced or the like) in the MRI apparatus, and the 
MR image is displayed. 

At Step Q2, an operator specifies a two-dimensional region of interest Rl 
(or a three-dimensional volumetric region of interest) on a displayed MR image 
Gl, as shown in Figure .4. 

At Step Q3, regular grid points gl, g2, g3, . . . are generated in the region of 
interest Rl (or in the volumetric region of interest), as shown in Figure 5. 

At Step Q4, points obtained by randomly moving the grid points gl, g2, g3, 
... in a two-dimensional (or three-dimensional) manner are defined as tracking 
start points SI, S2, S3, as shown in Figure 6. Random numbers for the 
random moving can be generated using a distribution function such as a 
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Gaussian distribution or uniform distribution. The flow then proceeds to Step 
Q5 in Figure 11. 

At Step Q5 in Figure 11, one of the tracking start points is selected. 

At Step Q6, diffusion tensor analysis is performed on the selected tracking 
5 start point in the three-dimensional image data collected by the diffusion tensor 
method in the MRI apparatus to determine the direction of the principal axis 
vector, the FA value, and the eigenvalues, and the principal axis vector is defined 
as a tracking direction vector. 

At Step Q7, if three-dimensional image data corresponding to a point at a 
10 unit distance along the direction of the tracking direction vector is present, the 
point is defined as a neighbor point and the flow proceeds to Step Q8; and if no 
three-dimensional image data corresponding to a point at a unit distance along 
the direction of the principal axis vector is present, the flow proceeds to Step Qll. 

At Step Q8, data at the neighbor point is created by interpolation or the like 
15 on the three-dimensional image data, and diffusion tensor analysis is performed 
to determine the eigenvectors, FA value, and eigenvalues. 

At Step Q9, if the FA value is equal to or more than a threshold, the flow 
proceeds to Step Q10 to continue the fiber tracking because the fiber tracking has 
not reached an end portion of a brain white matter fiber; and if the FA value is 
20 less than the threshold, the flow proceeds to Step Qll to terminate the fiber 
tracking because an end portion of a brain white matter fiber has been reached. 

At Step Q10, 

d i+1 = {XI (ei * di) ei + X2 (e 2 • di) e 2 + A3 (e 3 ■ di) e 3 } 
/ | XI (ei - di) ei + X2 (e 2 • di) e 2 + A3 (e 3 • di) e 3 1 
25 is set, where Al, X2 and A3 represent the eigenvalues of the diffusion tensor at a 
neighbor point, el, e2 and e3 represent the eigenvectors thereat, d*+i represents 
the tracking direction vector thereat, and di represents the tracking direction 
vector at the immediately preceding neighbor point or tracking start point. 

Figure 13 is a conceptual diagram showing the tracking direction vector 

30 di+i. 
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The flow then goes back to Step Q7. 

In this way, Steps Q7 — Q10 are repeated until no more three-dimensional 
image data are f ound or the fiber tracking has reached an end portion of a brain 
white matter fiber, and a fiber is tracked from the tracking start point SI to a 
5 neighbor point Nl, N2, N3, as exemplarily shown in Figure 7. At that time, 
connectivity is decided by using a scalar product of vectors, for example. 

At Step Qll, points from the tracking start point to the last neighbor point 
are saved as one brain white matter fiber. 

At Step Q12, if any tracking start point not selected at Step Q5 remains, the 
10 flow goes back to Step Q5; otherwise, proceeds to Step Q14 in Figure 12. 

At Step Q14 in Figure 12, an image of the saved brain white matter fibers as 
viewed in a desired view direction is produced, as exemplarily shown in Figure 
8. 

At Step Q15, the opacity at the tracking start point is defined as Xo. 
1 5 Moreover, 

X n +1 = FA n • X n 

is set, where X^+i represents the opacity at a neighbor point, FA n represents the 
FA value at the immediately preceding neighbor point or tracking start point, 
and X n represents the opacity thereat. 

20 At Step Q16, the display color (R, G, B) is defined as: 

R : G : B = 1 : A2 / XI .X3 / XI, 
where XI, XI and A3 represent the eigenvalues of the diffusion tensor. 

At Step Q17, the image of the fibers is displayed using the opacity X and 
the display color (R, G, B). 

25 According to the MRI apparatus of the second embodiment, the following 

effect can be obtained in addition to those in the first embodiment: 
(4) As shown in Figure 14, if the immediately preceding tracking direction 
vectors di and dj are different, the tracking direction vectors di+i and dj+i will be 
different even if the neighbor points Ni+i and Nj+i coincide with or He close to 

30 each other. Therefore, nerve fibers in different connection directions can be 
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distinguished based on the preceding connection directions even at a nerve fiber 
intersection C at which the nerve fibers in different connection directions 
intersect each other, as shown in Figure 15, and the nerve fibers can be tracked 
without mistaking the direction. 

In addition, to determine the tracking direction vector, an appropriate 
number N may be given to take an average vector of the next through N-th 
preceding tracking direction vectors.^ 
- Third Embodiment - 

Figure 16 is a flow chart showing fiber rendering processing by an MRI 
apparatus in accordance with a third embodiment. 

At Step Ql, an MR image in an axial or oblique plane is produced from 
three-dimensional image data collected by a diffusion tensor method or another 
imaging method (Tl- or T2-enhanced or the like) in the MRI apparatus, and the 
MR image is displayed. 

At Step Q2', an operator specifies a two-dimensional start region of interest 
Rl (or a three-dimensional start volumetric region of interest) and a 
two-dimensional end region of interest R2 (or a three-dimensional end 
volumetric region of interest) on a displayed MR image Gl, as shown in Figure 
19. 

At Step Q3, regular grid points gl, g2, g3, ... are generated in the start 
region of interest Rl (or in the start volumetric region of interest), as shown in 
Figure 5. 

At Step Q4, points obtained by randomly moving the grid points gl, g2, g3, 
... in a two-dimensional (or three-dimensional) manner are defined as tracking 
start points SI, S2, S3, as shown in Figure 6. Random numbers for the 
random moving can be generated using a distribution function such as a 
Gaussian distribution or uniform distribution. The flow then proceeds to Step 
Q5 in Figure 17. 

At Step Q5 in Figure 17, one of the tracking start points is selected. 

At Step Q6, diffusion tensor analysis is performed on the selected tracking 
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start point in the three-dimensional image data collected by the diffusion tensor 
method in the MRI apparatus to determine the direction of the principal axis 
vector, the FA value, and the eigenvalues, and the principal axis vector is defined 
as a tracking direction vector. 
5 At Step Q7, if a point at a unit distance along the direction of the tracking 

direction vector falls within the three-dimensional image data space, the point is 
defined as a neighbor point and the flow proceeds to Step Q8; and if the point 
falls outside the three-dimensional image data space, the flow proceeds to Step 
Qll. 

10 At Step Q8, data at the neighbor point is created by interpolation or the like 

on the three-dimensional image data, and diffusion tensor analysis is performed 
to determine the eigenvectors, FA value, and eigenvalues. 

At Step Q9, if the FA value is equal to or more than a threshold, the flow 
proceeds to Step Q10 to continue the fiber tracking because the fiber tracking has 
15 not reached an end portion of a brain white matter fiber; and if the FA value is 
less than the threshold, the flow proceeds to Step Qll to terminate the fiber 
tracking because an end portion of a brain white matter fiber has been reached. 
At Step Q10, 

d i+ i = {XI (ei - di) ei + XI (e 2 - di) e 2 + X3 (e 3 • di) e 3 } 
20 / | XI (ei • di) ei + XI (e 2 • di) e 2 + >3 (e 3 • di) e 3 | 

is set, where XI, XI and X5 represent the eigenvalues of the diffusion tensor at a 
neighbor point, el, e2 and e3 represent the eigenvectors thereat, d i+ i represents 
tracking direction vector thereat and di represents the tracking direction vector 
at the immediately preceding neighbor point or tracking start point. 
25 Figure 13 is a conceptual diagram showing the tracking direction vector 

di+i. 

The flow then goes back to Step Q7. 

In this way, Steps Q7 — Q10 are repeated until no more three-dimensional 
image data are found or the fiber tracking has reached an end portion of a brain 
30 white matter fiber, and a fiber is tracked from the tracking start point SI to a 
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neighbor point Nl, N2, N3, as exemplarily shown in Figure 7. At that time, 
connectivity is decided by using a scalar product of vectors, for example. 

At Step Qll, points from the tracking start point to the last neighbor point 
are saved as one brain white matter fiber. 

At Step Q12, if any tracking start point that not selected at Step Q5 remains, 
the flow goes back to Step Q5; otherwise, proceeds to Step Q13 in Figure 18. 

At Step Q13 in Figure 18, a decision is made as to whether the obtained 
fiber has an intersection with the end region of interest R2 (or the end volumetric 
region of interest), and the fiber is selected only if it has an intersection. 

At Step Q14', an image of only the selected brain white matter fibers / as 
viewed in a desired view direction is produced, as exemplarily shown in Figure 
20. 

At Step Q15, the opacity at the tracking start point is defined as Xo. 
Moreover, 

Xn+l = FA n * X n 

is set, where X n +i represents the opacity at a neighbor point, FA n represents the 
FA value at the immediately preceding neighbor point or tracking start point, 
and X n represents the opacity thereat. 

At Step Q16, the display color (R, G, B) is defined as: 
R : G : B = 1 : A.2 / XI : A3 / Al, 
where Al, A2 and A3 represent the eigenvalues of the diffusion tensor. 

At Step Q17, the image of the selected fibers is displayed using the opacity 
X and the display color (R, G, B). 

At Step Q18, a total sum with respect to all the selected fibers: 
MJValue = E Al • FA / L 
is calculated and displayed, where Al represents the first eigenvalue of the 
diffusion tensor of the selected fiber, FA represents the FA value thereof, and L 
represents the total length of the fiber. 

According to the MRI apparatus of the third embodiment, the following 
effects can be obtained in addition to those in the second embodiment: 
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(5) Since only the nerve fibers / passing through two sites are rendered, 
connectivity of the nerve fibers between the two sites can be visually recognized; 
and 

(6) Quantitative assessment is enabled by employing M_Value as an indicator 
5 of the strength of connection by nerve fibers between two sites. 

It is possible to display an average M_Value by dividing MJValue by the 
number of selected fibers. 

Moreover, the fibers may be displayed with the display brightness or 
display color changed according to MJValue. 
10 Many widely different embodiments of the invention may be configured 

without departing from the spirit and the scope of the present invention. It 
should be understood that the present invention is not limited to the specific 
embodiments described in the specification, except as defined in the appended 
claims. 
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